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Abstract. We study the implications of LHC searches on SUSY particle spectra using flat scans of the 
19-parameter pMSSM phase space. We apply constraints from flavour physics, — 2, dark matter and 
earlier LEP and Tevatron searches. The sensitivity of the LHC SUSY searches with jets, leptons and 
missing energy is assessed by reproducing with fast simulation the recent CMS analyses after validation on 
benchmark points. We present results in terms of the fraction of pMSSM points compatible with all the 
constraints which are excluded by the LHC searches with 1 fb~^ and 15 fb~^ as a function of the mass of 
strongly and weakly interacting SUSY particles. We also discuss the suppression of Higgs production cross 
sections for the MSSM points not excluded and contrast the region of parameter space tested by the LHC 
data with the constraints from dark matter direct detection experiments. 

PACS. 11.30.Pb Supersymmetry - 12.60.Jv Supersymmetric models 



1 Introduction 

Supersymmetry (SUSY) has emerged over the past two 
decades as possibly the best motivated model of new physics 
beyond the Standard Model (SM) . Together with stabilis- 
ing the masses at the electroweak scale, it provides gauge 
coupling unification and viable candidates for cold dark 
matter. Expectations for an observation of supersymmet- 
ric partners of the strongly interacting SM particles in the 
early stage of the LHC run have been high. In fact, global 
fits to constrained SUSY models, such as the CMSSM, 
including data from flavour physics, lower energy experi- 
ments and relic dark matter density have favoured SUSY 
particle masses below, or around, 1 TeV [IHl]- Recent 
searches for supersymmetric particles by the ATLAS and 
CMS experiments at the LHC have now excluded most 
of this portion of the parameter space for these mod- 
els, raising questions about the range of SUSY particle 
masses still allowed by the present data. Several stud- 
ies have been carried out to evaluate the impact of LHC 
and other data on SUSY parameters. Most of these stud- 
ies considered the highly constrained models with a small 
number of free parameters and large correlations between 
the masses of supersymmetric particles, which had been 
used for earlier benchmark studies and model parameter 
fits [3 [3]. In these constrained models, the LHC searches 
have resulted in a significant exclusion of masses at, or 



beyond, 1 TeV for the majority of SUSY particles. How- 
ever, these models are not representative of a generic min- 
imal supersymmetric extension of the Standard Model 
(MSSM) , where the particle mass parameters are indepen- 
dent. The phenomenological MSSM (pMSSM), with its 
19-parameter phase space [7j, was proposed to reduce the 
theoretical prejudices of these constraints. Two pioneer- 
ing studies considered the pMSSM for a global Bayesian 
fit using low energy and dark matter constraints to make 
inferences on SUSY particle masses j8| and the sensitivity 
of LHC searches [5^11j. More recently, two further stud- 
ies addressed the sensitivity of LHC to pMSSM scenarios 
based on ATLAS [H] and CMS [T3] analysis strategies. 

In this paper we report a study of the impact of the 
present LHC data at 7 TeV in the pMSSM with neutralino 
dark matter. Compared to the previous pMSSM study of 
Ref. |TT], we use here a larger scan statistics of almost 25M 
points associated to a broad range of the parameters, in 
particular for the gluino mass, in a flat scan, we adopt 
the selection criteria and results of the latest analyses on 
1 fb~^ of LHC data and the updated constraints from rare 
B decays, — 2 and ^?dm^^ and present our results in 
terms of fractions of pMSSM points excluded by the LHC 
data for 1 fb~^ and a 15 fb^-'^ projection without adopting 
a likelihood weighting. The main purposes of this study 
are to find which information on non-constrained SUSY 



scenarios are iraplied by LHC searches, and to investigate 
whether there are parameter regions not probed within 
the broad bounds on squark and gluino masses, possibly 
corresponding to distinct mass or couphng patterns. In 
addition, we aim to understand to which extent the spec- 
trum of weakly interacting SUSY particles, which are not 
directly probed by the LHC searches, is pushed towards 
higher masses as a consequence of the LHC results. 

This paper presents the scan and analysis method with 
results in the most studied neutralino dark matter sce- 
nario. To achieve our objectives, we perform flat scans 
over the parameter space of the pMSSM in order to probe 
different regions, without being interested in determining 
the best fit points of the model. Within these regions, we 
contrast the LHC sensitivity with results on dark mat- 
ter direct detection and we comment on the suppression 
of event yields in gg ^ hP ^ 77 channel which is most 
relevant for the search for light Higgs boson at the LHC. 

This work will be followed by analyses of the scan data 
in the framework of other scenarios, including light neu- 
tralino dark matter, gravitino dark matter, the NMSSM 
and in relation to the Higgs sector. 

This paper is organised as follows. Section [2] presents 
the scan of the SUSY parameters discussing the techni- 
cal tools and the constraints from lower energy data and 
dark matter relic density adopted for our analysis. The 
constraints from the LHC data are implemented by gen- 
erating event samples for each SUSY point and applying 
the analyses performed by the CMS collaboration to these 
events after fast simulation. Section |3] discusses the imple- 
mentation of these analyses in our scan framework and 
their validation using the official CMS simulation and re- 
construction. In section |4] we discuss the scan results in 
terms of both the sensitivity of the LHC searches with 
1 fb~^ and 15 fb~^ to strongly- and weakly-interacting 
sparticles. These results are relevant to the assessment of 
the physics potential of the 7 TeV LHC operation and 
also to the planning of a future lepton collider, which 
should operate at a centre-of-mass energy above the gaug- 
ino and slepton pair production threshold to perform de- 
tailed studies of the spectroscopy of the weakly interacting 
SUSY particles, if SUSY is indeed realised in nature. Fi- 
nally we discuss the spin independent XP scattering cross 
section relevant to dark matter direct detection experi- 
ments, the suppression of light Higgs production and de- 
cay yields compared to the Standard Model and the fine 
tuning for the points not excluded by the LHC data. 



2 pMSSM Scans 

In the pMSSM, all the soft SUSY breaking parameters 
are assumed to be real, in order to avoid new source of 
CP violation. The sfermion mass matrices and the tri- 
linear coupling matrices are taken to be diagonal, to sup- 
press FCNCs at tree level. In addition, the first and second 
sfermion generations are universal at low energy and the 
trilinear couplings are set to be equal for the three gener- 
ations. With these assumptions, we are left with 19 free 
input parameters at the weak scale. 



Parameter 


Range 


tan P 


[1, 60] 


Ma 


[50, 2000] 


Ml 


[-2500, 2500] 


M2 


[-2500, 2500] 


M3 


[50, 2500] 


Ad = As ^ Ab 


[-2000, 2000] 


~- A-c ^ A-f 


[-2000, 2000] 


A-e ^ 


[-2000, 2000] 


fl 


[-1000, 2000] 


Mer^ = Mfi^ 


[50, 2500] 


Me-^ = M^^ 


[50, 2500] 




[50, 2500] 




[50, 2500] 




[50, 2500] 




[50, 2500] 


Maj, = Mg^ 


[50, 2500] 




[50, 2500] 


= M~s, 


[50, 2500] 




[50, 2500] 


SUSY parameter ranges (in GeV when ap 


Parameter 


Value 


as(Mz) 


0.1184 


fhb{fhb) 


4.19 GeV 


pole 
mf 


172.9 GeV 



Table 2. SM parameter values [Ti]. 



To explore the pMSSM parameter space, we perform 
flat scans by randomly varying the 19 parameters within 
the intervals given in Table [T] with the SM parameters of 
Table [2j Our scan strategy, which is not guided through 
the parameter space as in the case of Markov Chain Monte 
Carlo methods, may miss local features. This does not 
change in a significant way the fraction of excluded points 
and the other results given in Section |4] as long as these 
regions are narrow. On the other hand our approach en- 
sures a rather uniform sampling of a very broad parameter 
space, without the risk of missing sizeable but largely dis- 
connected acceptable regions of parameters. In section |4] 
we discuss our scan coverage in terms of SUSY spectral 
features and particle properties. We also study how the 
fraction of valid MSSM points excluded by the LHC anal- 
ysis changes with the SUSY parameter ranges adopted for 
the scan. We generate a total of 24.57 million points. In 
this study, it is imposed that the lightest supersymmetric 
particle is the neutralino. 

For the flavour related quantities, we adopt the numer- 
ical inputs of [15J. 



2.1 Simulation Tools 

The study of the LHC sensitivity in the pMSSM param- 
eter space requires the use of several programs from the 
spectrum generation to the computation of flavour and 
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other low energy and dark matter related observables, as 
well as the simulation and analysis of events in 7 TeV pp 
collisions. Throughout the analysis we use the SUSY Les 
Houches Accord (SLHA) format [Hj to store the SUSY 
parameters, physical spectra and couplings, the HEPEVT 
format for simulated events and the ROOT format [T7| for 
results of the scans, the physics objects after fast simula- 
tion and the analysis results. 

SUSY spectra are generated with SOFTSUSY 3.1.7 [T8]. 
Results have been cross checked with SUSpect [7| and 
found to be compatible with the exception of few regions of 
the parameter phase space. The masses, widths and decay 
channels of Higgs bosons are obtained using FeynHiggs 
2.8.4 [13 [50] and HDECAY 3.53 [H]. The widths and 
branching fractions of the SUSY particles are computed 
using SDECAY 1 . 3 [22 . Flavour observables, the muon anoma- 
lous magnetic moment and the dark matter relic density 
are calculated with Superlso Relic v3.1 [TSllS^. Also, 
the neutralino-nucleon scattering cross-sections which are 
needed to evaluate dark matter direct detection constraint 
are computed with micrOMEGAs 2.4 [21]. Finally, we con- 
front the spectrum of each pMSSM point to the exclusion 
bounds obtained from Higgs at LEP, Tevatron and LHC 
using HiggsBounds 3.4.0 [^[Mj . 

We use the Superlso program as the central control 
program, which provides interfaces with the other above- 
mentioned programs through SLHA files. For each point 
giving a valid spectrum with Xi LSP, we store the pMSSM 
parameters, SUSY particle spectra, branching fractions, 
and flavour and low energy observables in ROOT format. 
We then impose the constraints on these observables and 
on the masses of the supersymmetric particles imposing 
the exclusion bounds from SUSY searches at LEP and 
the Tevatron and generate an SLHA file for each accepted 
point. These SLHA files, extended to include SUSY branch- 
ing fractions, are used as input to PYTHIA 6 .4. 24 [27] for 
event generation of inclusive SUSY production in pp inter- 
actions at 7 TeV, using CTEQ5L parton distribution func- 
tions [21] ■ Cross sections are rescaled to their NLO values 
with the k-factors computed for each point by Prospino 
2 . [2S] • Generated events are passed through fast detec- 
tor simulation using Delphes 1 . 9 [SC tuned for the CMS 
detector. The output of the Delphes simulation is saved in 
ROOT format and followed by a custom implementation of 
the CMS SUSY analyses. Results are stored in a database 
which is subsequently read-in in ROOT and used to deter- 
mine the observability of each of the pMSSM scan points. 
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Fig. 1. Branching fractions of _B — >■ Xs'y (left panel) and Ba — > 
/i^^^ (right panel) of pMSSM points as a function of the gluino 
and A'^ boson masses, respectively. The white horizontal lines 
correspond to the SM values. 



the experimental accuracy is good thanks to the B factory 
data. This branching fraction provides strong constraints 
on the SUSY parameter space, especially for large tan^, 
where it receives large enhancements. We use the following 
95% C.L. interval [ftl: 



2.16 X 10^'' < BR(B ^ X,j) < 4.93 x 10" 



(1) 



The decay Bg — > A'^M" is also a loop mediated process, 
which can receive extremely large SUSY contributions at 
large tan/3, and its branching fraction can be enhanced by 
several orders of magnitude compared to its SM prediction 
[57WT] . We use the two loop calculations based on [i^HiS] . 
This decay mode has not yet been observed, but recently 
LHCb and CMS provided upper limits on its branching 
fraction which are only 3.4 times above the SM value |46) . 
The CMS and LHCb combination yields a 95% C.L. upper 
limit of 1.1 X 10~*. Here we add theoretical uncertainties 
and use a constraint of [1T1: 



2.2 Constraints 



BR{Bs A^^M") < 1.26 X 10' 



(2) 



We apply constraints from flavour physics, the anoma- 
lous muon magnetic moment, relic dark matter and SUSY 
searches at LEP and the Tevatron. First we consider the 
flavour observables. The branching ratio of B ^ Xg'y has 
been thoroughly studied for its high sensitivity to new 
physics effects (see for example [3Tti34] ) , since the SM con- 
tributions only appear at loop level. We follow the NNLO 
calculation of [211 [33] for the theoretical predictions. The 
theoretical uncertainties here are well under control and 



The range of branching fractions of these processes in 
the pMSSM parameter space is displayed as a function of 
the gluino and the A'^ boson mass in Figure [ij 

In addition to these rare processes, we consider a set 
of tree-level observables, which are very sensitive to the 
charged Higgs mass as well as to tan (3. These include 
BR{B ^- Tiy), BR{D, tv), BR{B Dtv) and R^23 
which is related to the K — > jiv decay. We use the fol- 
lowing 95% C.L. intervals, which include theoretical and 
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experimental uncertainties [HJ |47l |48] : 



0.56 < 



BR{B Tv) 



< 2.70 



4.7 X 10^2 < BR(L>, ^ Tiy) < 6.1 x 10"^ ^ (4) 
2.9 X 10"^ < BR(B ^ D°Ti^) < 14.2 x 10"^ , (5) 
0.985 < Rf,23{K fii^) < 1.013 . (6) 

We should emphasise here that the allowed intervals 
for the flavour observables are quite robust within the 
MSSM and we have taken into account consistently the 
theoretical and experimental uncertainties. 

The anomalous magnetic moment of the muon, — 2, 
represents an outstanding open issue in precision low- 
energy measurements sensitive to SUSY contributions. The 
discrepancy between its SM prediction and the experimen- 
tal value obtained by E821 [JSIIS^, can be interpreted as 
a result of SUSY contributions [51], which would favour 
light masses of scalar muons ruling out large region of 
the SUSY parameter space [HllHlj. However, recent re- 
analyses of the hadronic contributions raise questions of 
the interpretation of the apparent discrepancy and 
offer a mean to reconcile the experimental result with SM 
prediction. Given the remaining uncertainty around the 
interpretation of the result, we adopt here the broad- 
est possible interval, which takes separately the lowest and 
highest 95% C.L. bounds given in the recent literature: 



- 2.4 X 10"^ < (5a^ < 4.5 x 10' 



(7) 



The next observable we use as a constraint is the dark 
matter relic density i7DM^^, for which we adopt the lat- 
est WMAP result [S7] from the cosmic microwave back- 
ground analysis. It has been shown [SSI HH] that many 
cosmological phenomena can strongly alter the value of 
the calculated relic density, leading to a very broad inter- 
val, 10"''^ < J7dm^^ < 10^, which takes into account un- 
certainties from realistic cosmological scenarios ^60|. Nev- 
ertheless, for this study we adopt the standard model of 
cosmology, and we consider that the lightest neutralino 
can account for the observed dark matter density from 
at least 0.1% of the WMAP value, up to the WMAP 
central value increased by 3.5 times the WMAP uncer- 
tainty summed in quadrature with 10% additional contri- 
bution to include theoretical uncertainties, in particular 
from the QCD equations of state in the relic density cal- 
culation [BT] . We therefore retain the following interval in 
our numerical analysis: 



10-'' < f2^h'^ < 0.155 



(8) 



In order to test the sensitivity of our analysis to the ad- 
doptcd range of these constraints, we increased the range 
of Eqs. ([l])-([6]) by a factor 1.5 and repeat the analysis. 
Results are discussed in Sections |4.1| and |4.2| In addition, 
the dark matter direct detection scattering cross-sections 
also provide relevant information. The CDMS [62 and 
the Xenon 100 collaborations have recently published 
stringent limits on the spin-independent elastic WIMP- 
nucleon scattering cross section. In this study, we calculate 



CMS Expected UL{1.1 fb^) 
This Study Expected UL (1.1 fb ^) 
Ttiis Study Expected UL (15 fb ^) 
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Fig. 2. 95% C.L. upper limit exclusion contours in the 
(mA,tan/3) plane. The dashed line shows the expected exclu- 
sion from CMS with 1.1 fb"'^ [64], the solid line the correspond- 
ing contour obtained with our simulations and the dotted line 
the same contour with 15 fb^^. 



the scattering cross sections for each SUSY model point, 
but we do not apply it as a constraint. We discuss the im- 
plications of dark matter direct detection on our results 
in Section [331 

Furthermore, we consider the constraints from direct 
Higgs searches at LEP, Tevatron and LHC, using the code 
HiggsBounds, which incorporates also the recent results 
on light and heavy Higgs searches at LHC. In order to 
study the scaling of heavy Higgs exclusion in {Ma, tan f3) 
plane we estimate the sensitivity of the CMS analysis in 
the TT channel [M]. We calculate the product ct^BRt-t- as 
a function of Ma and tan /3 in the NUHM model where we 
also vary mp, ^1/2, ^0 s-nd We reconstruct the wedge 
for the 95% upper bound for 1.1 fb^^ on tan/3 as a func- 
tion of Ma based on the expected sensitivity of the CMS 
analysis as given in [M]. We then rescale the 95% upper 
bound assuming an integrated luminosity of 15 fb"^. The 
results are shown in Figure [2] We use this c onst rain t for 
our study with 15 fb^^ of data in Sections 4.3 and 4.4 



Finally we impose the constraints on the SUSY masses 
summarised in Table |3] In the present study we consider 
for simplicity a lower limit of 46 GeV for the lightest neu- 
tralino mass, even if this limit can be circumvented in 
specific cases. 

The effect of these constraints on the masses of some 
SUSY particles is shown in Figures [3] and [4] where the 
fraction of valid MSSM points fulfilling the subsequent 
cuts on fiavour, — 2 and Q^h^ observables is shown as a 
function of the masses of g, hghtest q, xf and £ {i — e, fj,). 
We observe that the flavour cuts disfavour masses below 
~500 GeV for the gluino and the lightest scalar quark due 
to the deviations arising from their contributions to heavy 
fiavour loop-mediated decays. Very light, as well as heavy 
charginos are equally disfavoured by the relic dark matter 
constraint, which slightly favours light sleptons. However, 
with the exception of gluinos and lightest q below 250 GeV 
and lightest charginos of more than 1 TeV, the acceptance 
of all these cuts appear quite flat with the mass of both 
strongly and weakly interacting SUSY particles. 
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Fig. 3. Fraction of valid MSSM points selected after tlie flavour (dotted line), — 2 (dashed line) and Sl^h^ (solid line) 
subsequent constraints as a function of the masses of strongly-interacting SUSY particles: the gluino (upper panel) and the 
lightest scalar quark (lower panel). 
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Fig. 4. Fraction of valid MSSM points selected after the flavour (dotted line), — 2 (dashed line) and Sl^h^ (solid line) 
subsequent constraints as a function of the masses of two weakly-interacting SUSY particles: the lightest chargino xt (upper 
panel) and the scalar lepton (lower panel). 



3 LHC analysis simulation 

The sensitivity of the LHC SUSY searches with jets, lep- 
tons and missing energy for the pMSSM points selected 
from our scans has been assessed by reproducing on fast 
simulation the event selection criteria of three analyses 
of the CMS Collaboration. These analyses are the SUSY 
searches in all hadronic events with ut |65| . in same-sign 
isolated dilepton events with jets and missing energy [SB] 
and in opposite-sign dilepton events with missing trans- 
verse energy [SJ. Since these analyses are cut based, the 
observation or exclusion of a given model is obtained from 
the number of events observed in the signal region and 
the events selected by each analysis and signal regions are 
fully uncorrelated, they are well suited for a fast simula- 
tion and reconstruction study. While we have adopted the 
CMS analyses for this study, the ATLAS experiment has 
comparable sensitivity to SUSY. 

Events generated with PYTHIA are subsequently recon- 
structed and classified following the procedure of the three 
CMS analyses from the physics objects obtained from the 
Delphes fast simulation. Jets are reconstructed using the 



anti-kt algorithm [68] . implemented in the Fast Jet pack- 
age (SS] ■ The hadronic ax analysis considers events fulfill- 
ing the selection requirements of ref. |65j binned according 
to their Ht value. The four bins: 

• 575 <Ht < 675 GeV, 

• 675 <Ht< 775 GeV, 

• 775 <Ht < 875 GeV, 

• Ht> 875 GeV, 

are used to search for the SUSY signal. The same-sign 
dilepton analysis considers high-pf leptons in the four search 
regions defined by 

• Ht> 400 GeV and MET > 120 GeV, 

• Ht> 400 GeV and MET > 50 GeV, 

• Ht> 200 GeV and MET > 120 GeV, 

• Ht> 80 GeV and MET > 100 GeV. 

The opposite-sign dilepton analysis looks to the high MET 
and high Ht search regions defined by 

• MET > 275 GeV and Ht > 300 GcV, 

• MET > 200 GeV and Ht > 600 GeV, 
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Conditions 


X"! 


46 




X2 


62.4 


tan /3 < 40 


1 

X3 


99.9 


tan /3 < 40 


-u 
Xi 


116 


tan /3 < 40 


v± 
Xi 




Lail n ^ ^U, lti~±: — llb^n ^ xj V 
Xi Xi 


Br 






Sl 


107 




fi 


81.9 


mf, — m^o > 15 GeV 


Ur 


100 


vTLuj, — rri-^o > 10 GcV 


ul 


100 


rriuj^ — m^o > 10 GeV 


tl 


95.7 


rnf — m-o > 10 GeV 

11 Xi 


O^R 




Hi, l'l'<-^0 ^ Vjrc V 
Xi 


dL 


100 


— m-0 > 10 GeV 

°-L Xi 




248 
220 


m^o < 70 GeV, mr - m^o > 30 GeV 
m^o < 80 GeV, m^^ - m^o > 30 GeV 


hi 


210 
200 
100 


m^o < 100 GeV, mr - m^o > 30 GeV 
m^o < 105 GeV, mr - m^o > 30 GeV 

Xi oi X.1 

mr — m-0 > 5 GeV 

"1 A-| 


9 


195 





Table 3. Constraints on the SUSY particle masses (in GeV) 
from searches at LEP and the Tevatron 1141. 



respectively. The expected number of background events 
from SM processes in each search region is taken from 
the published CMS results, obtained with full simulation 
and reconstruction, validated on data, and rescaled to the 
assumed integrated luminosity. The 95% confidence level 
exclusion of each SUSY point in presence of background 
only is determined using the CLs method |70]. The num- 
ber of events in each of the search regions of the CMS 
analyses described above is computed for the signal plus 
background and the background only hypotheses, where 
the signal represents the number of events in the search re- 
gions estimated for the SUSY point under test. The num- 
ber of events in each of the search regions is assumed to 
be uncorrelated. 



3.1 Simulation Validation 

In order to ensure that the sensitivity to signal events in 
our fast simulation and reconstruction is representative 
of that from the analysis of the real CMS data, a vali- 
dation of this procedure is performed comparing results 
of the fast simulation used in this analysis to those of 
signal events fully simulated and reconstructed with the 
CMS code. These events are chosen from the set of LM 
benchmark points, defined in the CMSSM model and used 
by the experiment for studies of the analysis sensitivities 
|71j . The fast simulation and reconstruction validation is 
performed in three stages. First, we compare the shape 
of event variables after event preselection. These are the 
number of reconstructed jets, the (Xt variable used to dis- 
criminate against QCD jets and the sum of the jet trans- 
verse energies Ht for the hadronic analysis; the pt of the 



two leading leptons, the dilepton invariant mass, the event 
missing transverse energy, MET, and variable for the 
dilepton channels. The shape of the distributions in the 
CMS simulation and reconstruction and in the Delphes 
fast simulation are in good agreement. Then, we compare 
the number of signal events selected in the different sig- 
nal regions for each LM point. In general, the agreement 
between CMS full simulation and reconstruction and the 
Delphes fast simulation is good, with differences in the 
observed signal event rates below 25%. Finally, we com- 
pare the exclusion contours obtained using our fast simu- 
lation in the (mo, TO1/2) plane of the CMSSM model for 
the three analyses with the expected exclusion regions re- 
ported by CMS. Again, the 95% C.L. exclusion contours 
for the mo and TO1/2 CMSSM parameters obtained with 
our simulation agree within 20% to those of CMS as shown 
in Figure [5j 

4 Results 

Results are presented for 1 fb~^ and 15 fb~^ of inte- 
grated luminosity, corresponding to the statistics reported 
by CMS to the EPS 2011 conference and that expected by 
the end of 2012, respectively. These are based on the anal- 
ysis of the 835k pMSSM points fulfilling all the selection 
criteria discussed above and for which events have been 
successfully simulated, selected out of 24.57M pMSSM 
generated points. The statistics and efficiency if the sub- 
sequent selection steps, from the generated points giving 
a valid SOFTSUSY spectrum through the mass, low-energy 
data and Q^h? constraints to those for which a sample of 
events for inclusive SUSY production at 7 TeV is success- 
fully simulated with PYTHIA are summarised in Table [4] 
We discuss results for strongly interacting SUSY particles 
{g and q) , to which the present LHC searches are directly 
sensitive, and for their weakly interacting companions {i 
and x), which are more indirectly probed at the LHC and 
are instead of major interest for the planning of a future 
high energy lepton collider. We present results in terms of 
the fraction of pMSSM poi nts c ompatible with the con- 
straints discussed in Section 2.2 (referred to as "accepted 
pMSSM points" in the following) which are excluded by 
the LHC searches as a function of the mass of the SUSY 
particle of interest. We study the scan coverage in terms of 
the characteristics of the SUSY spectra and particle prop- 
erties. We observe that in 33.2% of accepted points, the 
gluino is lighter than the lightest of the squarks of the first 
two generations. 13.3% of the points have the gluino heav- 
ier than any of the squarks. Compressed spectra are par- 
ticularly challenging for searches at the LHC. We observe 
that only 1.5% of our accepted points have Mg — M^o < 
50 GeV, while 3.7% have - M^o < 50 GeV, where q 
is lightest squark of the first two generations. Turning to 
the particle properties, the Higgs boson mass ranges from 
80 to 135 GeV. The lightest neutralino is bino-like in 6%, 
wino-like in 23%, and higgsino-like in 24% of the accepted 
points. 

Further, we study the sensitivity of the three analy- 
ses used. We observe that of the points excluded with 1 
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Fig. 5. 95% C.L. exclusion contours in tlie (mo, mi/2) plane for CMSSM points with tan/3 = 10, = and sign(/i) = +1 for 
1 fb^^ of data at 7 TeV. The contours obtained from our fast simulation analysis (black solid line) are compared to the bands 
of the expected exclusion limit from the CMS full simulation and reconstruction including uncertainties shown by the lighter 
coloured lines. 
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Efficiency 


Generated Points 


24.57M 


1 


1 


Valid Spectra 


9.41M 


0.383 


0.383 


LSP and 








Mass Limits 


2.62M 


0.278 


0.107 


Higgs Limits 


1.81M 


0.691 


0.074 


Flavour and — 2 


1.34M 


0.743 


0.055 




897k 


0.668 


0.037 


SuccessfuU 








Simulation 


835k 


0.931 


0.034 



Table 4. Scan statistics. Points for which a 7 TeV pp inclusive 
SUSY event sample is successfully simulated and reconstructed 
are used in the subsequent analysis. 



(15) fb-i of data, 96.2 (96.9)% are excluded by the fully 
hadronic analysis, 27.6 (30.0)% by the same sign lepton 
and 10.4 (16.1)% by the opposite sign lepton analysis. 3.8 
(3.1)% of the points are excluded by leptonic analyses but 
not by the hadronic analysis. 



4.1 Strongly Interacting Sparticle Spectra of Allowed 
pMSSM Points 

We consider the masses of the gluino, g and of the lightest 
scalar quark of the first two generations. Figure [6] shows 
the fractions of accepted pMSSM points, which can be 
excluded by the results of the CMS analyses on 1 fb~^ 
and 15 fb~^ of data. We observe that in the case of the 
gluino the LHC data can exclude more than 85% of the 
pMSSM points up to a mass of --520 GeV and 700 GeV 
for 1 fb^^ and 15 fb~^, respectively. The small fraction 
of points escaping exclusion below these mass limits have 
either small production cross section or small mass differ- 
ence between the gluino and the LSP, resulting in com- 
pressed SUSY spectra yielding a lower transverse energy 
for the final state, as already pointed out in [TT]. Above 
these masses, the sensitivity of the LHC data decreases 
rather sharply and disappears for Mg ^1.5 TeV. For the 
case of the lightest scalar quark of the first two genera- 
tions, the points which can be excluded by the LHC data 




2000 2500 
M(g) (GeV) 




1500 2000 2500 

M(lightest q) (GeV) 



Fig. 6. Fraction of valid MSSM points excluded at 95% C.L. 
by 1 fb~^ and 15 fb~^ of LHC data as a function of the masses 
of the gluino g (upper panel) and the lightest scalar quark of 
the first two generations (lower panel). 



account for more than 85% of the accepted pMSSM points 
up to mass values of 320 GeV and 510 GeV for 1 fb~^ and 
15 fb~^ of data, respectively, while the effect of the LHC 
data disappears above 1.0 TeV and 1.3 TeV. These curves 
provide us with generic mass limits for strongly interact- 
ing SUSY particles which extend the results obtained on 
specific models, such as the CMSSM. As a comparison, 
in the CMSSM the limit on the gluino mass for 1 fb~^ of 
7 TeV data is between 500 GeV and 1.1 TeV depending on 
mg. In general, there are set of pMSSM parameters cor- 
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Fig. 7. Fraction of valid MSSM points excluded at 95% C.L. 
by 1 fb^^ and 15 fb~^ of LHC data as a function of tlie inclusive 
SUSY cross section. 



responding to spectra which are not detectable with the 
CMS analyses discussed above. Figure[7]shows the fraction 
of vahd MSSM points excluded at 95% C.L. by 1 fb'i and 
15 fb~^ of LHC data as a function of the inclusive SUSY 
cross section for 7 TeV collisions. Even at values of the 
SUSY cross section in excess of 5 pb which accounts for 
19% of accepted points, approximately 20% of the pMSSM 
points are not observables because of the low transverse 
energy in the SUSY decay products. These points deserve 
a special attention in future LHC studies. pMSSM points 
with very compressed spectra having a mass difference be- 
tween the gluino or the lightest of the scalar quarks of the 
first two generations and the LSP below the jet pt cut 
used in the analyses account for 5% of the total accepted 
points and 6% of those not excluded but having a SUSY 
inclusive cross section larger than 5 pb. 

We assess the sensitivity of our results on the adopted 
range for the input SUSY parameters and the constraints 
by studying the change of the fractions of accepted pMSSM 
points excluded with 1 fb~^ of data and Mg < 500 and 
Mq < 400 GeV on a dedicated scan of IM generated 
points. These fractions are 97.8% and 89.9% for our stan- 
dard ranges and become 95.3% and 83.4% if we increase 
the SUSY parameter range by a factor 1.5. This shows that 
the results of these analyses in the range of sensitivity of 
the LHC data are only moderately sensitive to the chosen 
range of SUSY parameters. If we make looser cuts for con- 
straints ([l])-([6]) corresponding to 3.5 a range the fraction 
of points passing this cut increases from 76% to 85%, and 
the fractions above become 95.6% and 83.3% for gluinos 
and squarks respectively. Finally, we apply tighter relic 
dark matter density cut corresponding to the 95% C.L. of 
the seven year WMAP resuh (0.0924 < n^h'^ < 0.1316). 
This reduces the fraction of accepted points from 66.8% to 
0.7% and the fractions of excluded points in the selected 
gluino and squark mass range become 99.8% and 95.0%. 



4.2 Weakly Interacting Sparticle Spectra of Allowed 
pMSSM Points 

Likelihood analyses of constrained SUSY models have in- 
dicated that, as a result of the LHC searches, which have 
excluded the region of parameters favoured by data before 
the start of the LHC, the masses of SUSY particles over 
the allowed portion of the parameter space have shifted to 
larger values [72J . Instead, the impact of the LHC searches 
for gluinos and scalar quarks on weakly interacting su- 
persymmetric particles in the general MSSM is only indi- 
rect. In the absence of the mass relations between particles 
of the strongly- and weakly-interacting sectors, typical of 
constrained models such as the CMSSM or the NUHM, 
the masses of gauginos and sleptons in the MSSM are a 
priori uncorrelated to those of the gluino and the squarks, 
since their mass parameters are free and independent. 

Correlations are introduced either by the constraints 
applied from flavour physics and relic dark matter den- 
sity or by the signatures in the LHC analyses. A correla- 
tion between the xf and the q masses is observed in this 
study and originates through the dark matter relic density 
constraint. Correlations between the masses of the q and 
those of x^, x", arise from searches in final states with 
leptons and missing energy. These topologies arise from 
cascade decays involving gauginos and sleptons, such as 
Q QXi X — ^ -^^j — ^ ^Xi- The negative result in searches 
with these topologies suppresses pMSSM points with x 
and £ with large enough mass splittings from the q to give 
large transverse energies to the leptons searched in the 
analyses. As a result of these correlations, the negative 
result in the search for strongly interacting supersymmet- 
ric particles may allow us to make an inference on the 
masses of gauginos and sleptons in the allowed region of 
the parameter space. Results are again presented as the 
fractions of accepted pMSSM points, which can be ex- 
cluded by the results of the CMS analyses on 1 fb~^ and 
15 fb-i of data as a function of SUSY particle masses 
(see Figure |8|. As expected, the CMS analyses prefer- 
entially reduce the solutions with smaller masses for the 
sleptons and lighter gaugino states. However, ^-^50% of the 
accepted pMSSM points with masses below 400 GeV are 
not excluded within our scans. We also observe that the 
domain of SUSY weakly-interacting particle masses above 
500 GeV is virtually unaffected by the present LHC data. 
This result contrasts the indications obtained in highly 
constrained SUSY models. These observations are impor- 
tant both for motivating searches of direct gaugino pro- 
duction at the LHC and for interpreting the LHC data 
to guide the choice of the energy scale of a future lepton 
collider. 

Similarly to what done for strongly interacting spar- 
ticles, we assess the sensitivity of these results on the 
adopted range for the input SUSY parameters and the 
constraints. Here we study the change of the fractions of 
accepted pMSSM points, excluded with M^o < 400 GeV 
and Mj < 400 GeV on a dedicated scan. These fractions 
are 55.9% and 57.1% for our standard ranges with 1 fb^^ 
of data. They become 37.8% and 38.2% if we increase the 
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Fig. 8. Fraction of valid MSSM points excluded at 95% C.L. by 1 fb"^ and 15 fb"^ of LHC data as a function of the masses 
of Xi (upper left), Xi (upper right), X2 (mid left), X2 (mid right), X2 (lower left) and the lightest slepton of the first two 
generations £^ (lower right). 



SUSY parameter range by a factor 1.5, 38.1% and 37.9% 
if we loosen the constraints ([l])-([6]) to a 3.5 cr range, and 
66.4% and 68.9% for a relic dark matter density cut corre- 
sponding to the 95% C.L. of the seven year WMAP result. 
The large change with the range of SUSY parameters is 
expected since the LHC data has little or no sensitivity on 
the mass of these particles and the fractions change with 
the broadening of the parameter space. 



4.3 Dark matter direct detection 

The broad coverage of the pMSSM phase space through 
the scans provides us with an opportunity to contrast the 



impact of the LHC SUSY searches and the dark mat- 
ter direct detection experiments. In particular, the recent 
results by the CDMS P] and Xenon 100 P] collabo- 
rations have moved the sensitivity of these experiments 
well within the region characteristics of the MSSM with 
Xi LSP [75] . Figure |9] shows the distribution of the spin 
independent XP scattering cross section as a function of 
the LSP mass for all the accepted pMSSM points and for 
those points not excluded by 1 and 15 fb^^ of LHC data 
together with the 90% C.L. exclusion bounds from CDMS 
and Xenon. The current 90% C.L. exclusion contour of the 
Xenon 100 experiment cuts through this region with about 
20% of the accepted pMSSM points giving a cross section 
exceeding the Xenon 100 limit. We note that 1 fb~^ of 



9 




CDMS Soudan 
Xenon 100 

Xenon 2012 Projection 



St 1,°,:. ° 



200 400 600 SOO 1000 1200 1400 

H(x;) (GeV/c^) 



Fig. 9. xp scattering cross section as a function of the LSP mass. The black dots represent accepted pMSSM points, tlic grey 
dots the subset of points not excluded by 1 fb~^ and the light grey dots those not excluded by 15 fb^^ of LHC data. The lines 
represent the 90% C.L. exclusion contours set by the CDMS (dotted) and Xenon 100 (thick continuous) experiments together 
with the Xenon projection for the 2012 run (thin continuous). 
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Fig. 10. Fraction of accepted MSSM points not excluded by 1 fb~^ of LHC data but excluded at 90% C.L. by the Xenon 100 
direct detection data as a function of the masses of g (upper left), Xi (upper right), xt (lower left) and the lightest slepton of 
the first two generations 1^ (lower right). 



LHC data preferentially removes points with light Xi in 
the large xp scattering cross section region which is also 
incompatible with direct detection experiment data. Since 
the main contribution to the XP cross section comes from 



diagrams mediated by the boson, the neutralino dark 
matter direct detection exclusion region preferentially re- 
moves points at low Ma and large tan /3 values. However, 
this constraint does not significantly impact the mass dis- 



10 



tributions for other SUSY particles. The fraction of ac- 
cepted pMSSM points incompatible with the Xenon 100 
data is shown as a function of the mass of the gluino, 
lightest slepton of the first two generations, Xi and X2 in 
Figure [T0| 



4.4 Higgs couplings 



The search for the Higgs boson is one of the most im- 
portant themes of the present LHC physics program. The 
current data has narrowed the allowed region for a SM- 
like Higgs boson down to mass interval 114 GeV < Mh < 
145 GeV. This is a tantalising indication that the likely 
mass region for the Higgs boson exactly coincides with 
that favoured by electroweak precision data for a SM Higgs 
boson and by SUSY for the lightest Higgs state, The 
points in our scan which are not excluded by 1 fb^^ of 
LHC data have Mh in the interval 98 to 130 GeV. The 
average value of Mh for these points is 118.2 GeV and we 
do not observe any significant shift compared to the value 
for all the accepted points. By the end of the 2012 the 
LHC should deliver enough luminosity, for both ATLAS 
and CMS to exclude the existence of a SM-like Higgs bo- 
son, if no signal is observed. It is important to understand 
how this prediction translates for the case of h°. Its cou- 
plings to fermions and gauge bosons are shifted by SUSY 
corrections [74]. As a result, in SUSY both the production 
cross section gg — > hP and the branching fractions for the 
decay processes /i" — >■ 77 and W^W~ differ from their 
SM values. We quantitatively study these differences for 
the accepted pMSSM points which are not excluded by 
the CMS searches on 1 fb^^ and 15 fb^^ of statistics, as 
discussed above. The single most sensitive variable is the 
mass of the boson. 



Figure 11 shows the range of values taken by the ratio 



BR(/i°->77) 



>77) 



as a function of Ma for 



accepted pMSSM points and not excluded by the CMS 
searches for a heavy Higgs boson and for SUSY particles 
with 1 fb^^ and 15 fb~^ of data as a function of Ma- 
The LHC searches for strongly interacting particles do 
not significantly change the prediction of the SUSY sup- 
pression for the product of production cross section and 
decay branching fraction, i?. Figure [12] shows the distri- 
bution of R for all the accepted pMSSM points, for the 
points not excluded by the LHC searches for gluinos and 
scalar quarks and also for the points compatible with the 
Xenon 100 direct detection experiment. The 2011 -|- 2012 
statistics will enable tighter limits on the combination of 
Ma and tan /?, if no signal is observed. However, it is in- 
teresting to observe that pMSSM points with suppression 
of the product of production cross section and branching 
fraction compared to the SM value of a factor of two or 
larger will remain. Even in the decoupling regime, i.e. for 
Ma > 500 GeV and beyond the A° boson sensitivity of 
the LHC 7 TeV run, there are acceptable pMSSM solu- 
tions where R is as small as 0.75. 
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Fig. 11. Distribution of the ratio R — 

BR(h°^77) 



J" — T 



as a function of Ma for accepted pMSSM points 



BR(ffsM^77) 

not excluded by 1 fb~^ (dark grey) and 15 fb~^ (light grey) of 
LHC data. The shaded area shows the full range of the values 
of R over the accepted pMSSM points. 
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12. Distribution of the ratio R 



X 



B'R{H° "H77) ^'^^ accepted pMSSM points (solid line), the same 
and not excluded by 1 fb^^ of LHC data (dashed line), the 
same and not excluded by Xenon 100 data (dotted line) and 
the same and not excluded by 15 fb~^ of LHC data (dashed- 
dotted line). 



4.5 Fine tuning 



A quantitative way of measuring how natural a model 
point is, can be achieved by introducing the fine tuning 
parameter. We use here an analogous to the fine tuning 
measure introduced first in [75] and [76] , although we fol- 
low mostly the approach presented in [77] and focus on 
the amount of fine tuning at the weak scale, rather than 
at the GUT scale, of the MSSM parameters. 

The mass of the Z boson in the MSSM is given by 
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at tree level, where 
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cos 2/3 



(9) 



(10) 



Following [751 ES] I we quantify fine tuning by computing 



d log Ml 
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(11) 



where ^ — Mjj^^ , M|f^ , b and /j, are the relevant parameters. 
This gives: 
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The overall fine tuning is obtained by adding the four 
5's in quadrature: 



A^[S{f,r + S{bf + S{Mlf + SiM'^J 



21 1/2 



(12) 



The larger this value, the more fine tuned is the SUSY 
parameter point. As a reference, the well known SPSla 
point [75] (now excluded by LHC data) has a fine tuning 
of about 70 in this definition. Figure [13] shows the distri- 
bution of fine tuning of our accepted pMSSM points, of 
those not excluded by the LHC data and by the Xenon 100 
scattering cross section upper bound. The shape of the fine 
tuning distribution is modified by these constraints, which 
preferentially removes particles with low masses, shifting 
the distribution towards larger fine tuning values. How- 
ever, a significant population of pMSSM solutions at low 
to moderate values of fine tuning survives after these cuts. 



5 Conclusions 

The implications of LHC searches on SUSY particle spec- 
tra have been studied using fiat scans of the 19-parameters 
pMSSM phase space. We evaluate the sensitivity of flavour 
physics observables to SUSY particle masses, and apply 




2 2.5 

log,„(A) 



Fig. 13. Distributions of the fine tuning variable A after the 
flavour, 3^ — 2 and nuuh^ (solid line), LHC (dashed line) and 
Xenon 100 (dotted line) subsequent constraints. The last bin 
includes overflows. 



them as constraints together with (7^ — 2, dark matter and 
the mass bounds from earlier LEP and Tevatron searches. 
The sensitivity of the LHC SUSY searches with jets, lep- 
tons and missing energy is assessed by reproducing the 
recent CMS analyses with fast simulation, validated on 
benchmark points. We study the fraction of pMSSM points 
compatible with all the constraints which are excluded by 
the LHC searches with 1 and 15 fb~^ as a function of the 
mass of strongly and weakly interacting SUSY particles on 
a sample of 835k accepted points, without performing an 
importance weighting of the points. We observe that the 
LHC data can exclude more than 85% of these pMSSM 
points up to a gluino mass of ~520 GeV and 700 GeV for 
1 fb~^ and 15 fb~^, respectively. On the contrary, the mass 
spectra of most of colour singlet states are only weakly im- 
pacted and, for the range of SUSY parameters adopted in 
this study, ^50% of the accepted points with masses below 
400 GeV are not excluded. In the pMSSM, the domain of 
SUSY weakly-interacting particle masses above 500 GeV 
is virtually unaffected by the present LHC data contrary 
to the results in highly constrained scenarios. Compared 
to the study of [13J, which is based on a similar method- 
ology for assessing the SUSY sensitivity of the LHC with 
CMS analyses, the results reported here for the impact of 
1 fb^^ of LHC data are in agreement. 

Comparing to dark matter direct detection experiments, 
the 95% C.L. exclusion contour of the Xenon 100 exper- 
iment cuts through the region of the accepted pMSSM 
points with about 20% of them exceeding the Xenon 100 
bound. The LHC data preferentially removes points in 
the large XP scattering cross section. A significant fraction 
of pMSSM solutions compatible with both the LHC and 
dark matter direct detection experiment data have low 
to moderate values of the fine tuning parameter. Finally, 
we use our sample of accepted pMSSM points to estimate 
the suppression of the light Higgs production cross section 
and 77 decay branching fractions compared to the Stan- 
dard Model predictions. We find that a suppression of up 
to more than a factor of two is compatible both with the 
current LHC data and with that to be collected by the 
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end of 2012, which has imphcations on the Higgs boson 
detectabihty. 
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